Charge and spin separation in the ID Hubbard model 
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We consider the repulsive Hubbard model in one dimension and show the different mechanisms 
present in the charge and spin separation phenomena for an electron, at half filling and bellow half 
filling. We also comment recent experimental results. 
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Electrons possess charge and spin, as strong as this 
unified picture may be imprinted in our minds these two 
quantities have been discovered separately. It is possible 
to trace back the discovery of the electron charge from 
the first experiments of J.J. Thomson to the measure- 
ments of R. Millikan that demonstrated the existence of 
a fundamental electric charge. Only sixteen years later 
G. E. Uhlenbeck and S. A. Goudsmith discovered that 
the electron not only has a charge, but also carries spin. 
The electron was the first elementary particle discovered. 
It is hard to recall any discovery, since then, with such 
an impact, not only on physics but on science and tech- 
nology as well. 

Since the first proposals ^-|), of a model of a lo- 
cal repulsive interaction as an appropriate system for 
modelling magnetism and metal-insulator transitions the 
so called Hubbard model became the prototype for the 
study of strongly correlated electrons. 

Strongly correlated electronic systems are under in- 
tense investigation, due to their possible role in high- 
T c superconductivity. In particular, the two dimensional 
Hubbard model is in the very center of those investi- 
gations, but one cannot find an exact solution and many 
results must rely on numerical calculations. On the other 
hand the one dimensional Hubbard model is an integrable 
model and it is believed that some of its most impor- 
tant features carry over higher dimensions 0]. Also the 
most outstanding physical phenomenon happens in the 
ID Hubbard model. Charge separates from spin. 

The Hamiltonian of the one-dimensional Hubbard 
model on a periodic L-site lattice may be written as 
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and 



are canonical creation and annihilation 



Fermi operators, n,j ltT = c^ a Cj,a is the particle number 
operator for electrons of spin a in the Wannier state cen- 



tered around the site j and U is the coupling constant. 
Periodicity is guaranteed by setting cl+\,<t — ci iCr . The 
Hubbard Hamiltonian conserves the number of electrons 
A" and the number of down spins M, Lieb and Wu § 
solved it, for fixed A" and M by means of the nested 
Bethe ansatz ||. They calculated the ground state en- 
ergy at half filling, N/L = 1. A review of the main 
results in the Hubbard model can be found in [Q. Some 
of the important symmetries of this Hamiltonian are its 
invariance under particle-hole transformations, under re- 
versal of spins H, allowing one to set 2M < N < L. 
In the following we shall denote the positions and spins 
of the electrons by Xj and aj, respectively. The Bethe 
ansatz eigenfunctions of the Hubbard Hamiltonian (|l|) in 
a sector Q are given as 



ip(xi, . . . ,x N ;ai, . . . , crjv) 



Q is a permutation related to one of the AH possible 
orderings of the coordinates of N electrons. The P- 
summation extends over all permutations of the num- 
bers 1, ...,N, these permutations form the symmetric 
group Sn and the function sign(Q) is the sign function 
on the symmetric group, which is —1 for odd permuta- 
tions and +1 for even permutations. The spin dependent 
amplitudes tpp(<jQi, . ■ . , ctqjv) are given in Woynarovich 
B and have form the Bethe ansatz wave functions of an 
inhomogeneous XXX spin chain. Two sets of quantum 
numbers characterize the wave functions (^|) the charge 
momenta {kj} and the spin rapidities {A;}. They satisfy 
the Lieb-Wu equations 
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A/ — sin kj — if// 4 
A; - sin kj + if// '4 
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A; — sin kj — if//4 
A* - sin kj + iU/4 

..,N , 

Xi - A m - iU/2 
Xi - X m + iU/2 



(3) 
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The wave functions (B) are joint eigenfunctions of the 
Hubbard Hamiltonian (|lj) and the momentum operator 
with eigenvalues 



JV 



U 



3=1 
N 



-2^2coskj + —{L-2N) , 

p = fc j J mod 2n ■ ( 5 ) 

Another crucial property, for the study of the spec- 
trum, is the invariance of ([!]) under the so(4) symmetry 
||. The Hubbard Hamiltonian (|l|) is invariant under ro- 
tations in spin space. The corresponding su(2) Lie alge- 
bra is generated by the operators 

L L 
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L 

f^EKr^), (6) 
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where [C, C f ] = ~X* , ICC] = C and [C f ,C"] = "C f - 

When the length L of the lattice is even there is an- 
other representation of su(2), which commutes with the 
Hubbard Hamiltonian || and generates the so-called re- 
pairing symmetry, 

L L 

r, = ^2(-iyc jlCjl , ^ = ^(-l)^ c + c + 

3=1 3=1 
L 
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where [77, 77^] = — 2i] z , [rj,rj z ] = rj and [77^ , 77^] = — r/' . 

The role of the 77-pairing symmetry is to connect sec- 
tors of the Hilbert space with different numbers of elec- 
trons. The generators of both the rotations in spin space 
and 77-pairing symmetry commute with one-another and 
combine into a representation of su(2) ® su(2)/Z2 = 
so(4). Making use of this algebraic structure it was shown 
that the set of Bethe-ansatz eigenstates can be extended 
leading to a complete set of eigenstates jnj . 

As a result from the above so(4) symmetry it is possi- 
ble to classify the excitation spectrum over the half-filled 
ground state, with respect to its representations. The 
charge and spin degrees of freedom separate. The scat- 
tering states of only four elementary excitations provide, 
in the repulsive case, the whole spectrum of excitations. 
Two of these elementary excitations transform as a triplet 
under the spin transformations and as a singlet under the 
77-pairing. This means that they carry spin but no charge. 
The resulting representation can be denoted by (^, 0) and 



one calls them spinons. The (0, |) representation of so(4) 
is associated to the other two elementary excitations, the 
holon and the antiholon, which carry charge but no spin. 
Bellow half filling there are infinitely many excitations, 
although only few of them are essential for low energy 
physics. 

At half filled band the energy and momentum of the 
spinons are expressed in terms of a spectral parameter A, 
which plays a role similar to the rapidity for relativistic 
particles. This excitation is gapless and the dispersion 
curves for the energy e s (A) and momentum p s (A) are 



e,(A) = 2 



p s (A) = -tt/2 - 



dcj Ji(u>) cos(o;A) 
lo cosh(cj?7/4) 

du> Jo(u>) cos(wA) 
uj cosh(w[//4) 



(8) 



where < p s (A) < 7r and Jo,i are Bessel functions. 

The holon and antiholon momenta differ by 7r as 
does not commute with the momentum operator chang- 
ing it by it. Their dispersion curves |8| are parametrized 
by k h the bare momentum of the holon and the spectral 
parameter is now sin(k h ). For the holon, which has a 
gap, we have, 



U 



c/.(A) = 17 + 2cos(k h ) 



(9) 



+2 



Pfc(A) = tt/2 - k 



cosh(uU/4) 
h f°° dtu Jq(lu) cos(wA) 
uj cosh(o;C//4) 



Applying the operator the energy and momentum 
of the antiholon are found to be e a h(k) = £h(k),p a h{k) = 
-7T +p h (k). 

Together the Lieb-Wu equations and the so(4) symme- 
try give a complete description of all the excitations. The 
usual approach is to determine the Bethe ansatz excita- 
tion taking the zero-temperature limit of the thermody- 
namic Bethe ansatz, making use of the string hypothesis. 
The Bethe states are then parametrized by sets of spec- 
tral parameters {kj}, {A a }, {A a /}, which are subject to 
a set of coupled algebraic equations [plj . The scattering 
states of these excitations provide the two body scatter- 
ing matrix. In |l3f| it was explained how to deal with the 
non trivial structure of a ground state related to many- 
body wave functions of Bethe ansatz solvable models in 
order to obtain the scattering matrix. For the Hubbard 
model the complete scattering matrix is a block diag- 
onal 16 x 16 matrix, where each 4x4 block describes 
the different scattering states such as the, spinon-spinon, 
spinon-holon, holon-spinon and holon-holon. 

Let us now examine the charge and spin separation at 
the half filled band in more detail. We first remark the 
differences in the excitation spectrum in the free, U = 0, 



2 



and interacting, U > 0, cases. When U = the elec- 
tron is the only excitation, it carries charge and spin. 
In the interacting case the electron disappears from the 
spectrum and we have instead the spinons, which carry 
spin, and the holons, which carry charge together with 
the antiholons. When U > the low energy dispersion 
curve for the spinons is linear, indicating gapless exci- 
tations, but the holon and antiholon dispersion curve is 
quadratic an the excitation has a gap. There are no other 
excitations at any energy. After separation occurs the 
spinon and holon still interact. The scattering state of a 
spinon holon can be constructed by choosing N = L — 1 
and M = 4 — 1, which means that one hole is intro- 
duced in the kj's and A^'s distributions. The spin and 
repairing of this excitations are given by S = \ = r\ and 
S z = —\ = rf . The scattering states of two spinons and 
two holons can also be constructed Knowledge of 

those states allows the calculation of the scattering ma- 
trix S of a spinon with momentum p on the standing 
holon as, 



p-i2J Q {i$j) 



(10) 



In the above expression the term, 2Jo(*^j-), is the inter- 
action constant between the holon and spinon and Jo is 
a Bessel function. 

Bellow half filling the situation is quite diverse, one can 
calculate the spin and charge velocities and realize that 
they are different. Meaning that if we drop an electron 
on the chain the spin will fly away from the charge. From 
another point of view one can consider the spectrum of 
low lying excitations and verify that it is possible to clas- 
sify them into two categories, i) possessing charge but no 
spin and gapped or gapless, ii) possessing spin but no 
charge and gapless. The electron state being gapped. 

Yet another possibility to approach the charge and spin 
separation property is to consider correlation functions at 
finite temperature T: 



g{x,t,T) 



tr{e--c](x,t)ci(0,0)} 



(11) 



Long distance asymptotics of these correlation func- 
tions is governed by the conformal dimensions A^. In 
the case of the Hubbard model the conformal dimensions 
A^ s associated to the charge and spin distributions have 
been calculated in It is then possible to realize 

that, asymptotically, the correlation function factorizes 
into two factors a charge and a spin one, which we find 
to be given by, 



JJezp{- 



-2ttT 



Af\x =F v c t\}exp{- 



-2vrT 



A±|a:Tf fl t|}, 



where v c and v s are the Fermi velocities for the charge 
and spin excitations. This, in turn, implies the factor- 
ization of the creation operators of an electron c| a into 
spin and charge factors. 

We have thus discussed the phenomenon of charge and 
spin separation in the lZJHubbard model. With the 
present exact results we expect to provide a basis for 
the usual intuition developed in the studies of strongly 
correlated electron systems (see for instance fll5|). 

It is worth noticing that numerically charge and spin 
separation was shown in some extensions of the Hubbard 
model ]lq| . Also in the framework of the t — J model 
charge and spin separation occurs [ fl7| . In the two dimen- 
sional realm a theoretical proposal, exploring a discrete 
gauge symmetry, to investigate charge and spin separa- 
tion has been put forward by Senthil and Fisher [ jl8| . 

More recently, several groups started experimental 
verification of charge and spin separation using the 
technique of angle-resolved photoemission spectroscopy 
(ARPES) in quasi one-dimensional materials. ARPES 
experiments can provide detailed dispersion curves 
through the analysis of the energy and emission angle 
of electrons hit by monochromatic photons. Using this 
technique one group [|l9| studied the ID charge-transfer 
insulator SrCu02, another experimental group studied 
the ID Mott- Hubbard type insulator NaV 2 5 @. The 
results coming from both experiments require a theoreti- 
cal framework where the charge and spin degrees of free- 
dom separate. 
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